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Abstract 

Congenital long QT syndrome is a heritable family of arrhythmias caused by mutations in 13 genes encoding ion channel 
complex proteins. Mounting evidence has implicated the Purkinje fiber network in the genesis of ventricular arrhythmias. In 
this study, we explore the hypothesis that long QT mutations can demonstrate different phenotypes depending on the 
tissue type of expression. Using computational models of the human ventricular myocyte and the Purkinje fiber cell, the 
biophysical alteration in channel function in LQTl, LQT2, LQT3, and LQT7 are modeled. We identified that the plateau 
potential was important in LQTl and LQT2, in which mutation led to minimal action potential prolongation in Purkinje fiber 
cells. The phenotype of LQT3 mutation was dependent on the biophysical alteration induced as well as tissue type. The 
canonical AKPQ mutation causes severe action potential prolongation in both tissue types. For LQT3 mutation F1473C, 
characterized by shifted channel availability, a more severe phenotype was seen in Purkinje fiber cells with action potential 
prolongation and early afterdepolarizations. The LQT3 mutation S1904L demonstrated striking effects on action potential 
duration restitution and more severe action potential prolongation in Purkinje fiber cells at higher heart rates. Voltage 
clamp simulations highlight the mechanism of effect of these mutations in different tissue types, and impact of drug 
therapy is explored. We conclude that arrhythmia formation in long QT syndrome may depend not only on the basis of 
mutation and biophysical alteration, but also upon tissue of expression. The Purkinje fiber network may represent an 
important therapeutic target in the management of patients with heritable channelopathies. 
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Introduction 

Congenital long QT syndrome (LQTS) is a lieritable family of 
arrhythmias caused by mutations in 13 genes encoding ion 
channel complex proteins[l]. These mutations exert their effect 
througli altered ion channel gating and/ or trafficking, resulting in 
prolongation of action potential duration (APD). This prolonga- 
tion manifests as lengthening of the QT interval on the 
electrocardiogram, providing an electrophysiological substrate 
for arrhythmia formation. While the majority of the characteriza- 
tion of arrhythmias in this disorder to date has involved the 
ventricular myocyte (VM), mounting evidence implicates the 
cardiac specialized conduction system, comprised of the Purkinje 
fiber (PF) network, as a particularly important trigger of 
arrhythmias [2-6] . PF cells (PFC) have distinctive action potential 
(AP) characteristics, including the property of automaticity, longer 
APD, and lower plateau membrane potential [7]. Each of these 
features may uniquely predispose these cells to arrhythmia 
formation. 

Mathematical modeling has been employed to simulate genetic 
and acquired channelopathies and assay their efTects on cardiac 
tissue [8,9]. We have recentiy developed comprehensive integrative 



models of the electrophysiology of both the human VM[10] and 
PFC [11]. By virtue of their biophysicaUy detailed Markov 
representations of ion channels (constrained to data obtained 
from human ion channels in heterologous expression), these 
models are well-suited to the study of LQT mutations, through 
targeted modification of the processes responsible for channel 
dysfunction. 

In this study, we create models of LQTS in VM and PFC to 
investigate the hypothesis that LQTS-related mutations may 
manifest their phenotype in a tissue-specific manner. Several 
disease-causing mutations are simulated in each cardiac tissue 
type. We demonstrate that arrhythmia formation in LQTS may 
depend not only upon the specific mutation and biophysical 
alteration but also upon the tissue of expression, supporting an 
important role for the PFC. 

Results 

LQTl and LQT2 mutations have a mild phenotype In PFC 

As LQTl and LQT2 mutations often result in complete loss-of- 
function of the affected KCNQ_1 or HERG allele, respectively 
(through defective trafficking or biophysical changes in channel 
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function) [1], idealized heterozygous mutations were simulated by 
reducing peak conductance of IKs or IKr by 50% to simulate 
heterozygous LQTl or LQT2, respectively. Figure 1 shows the 
resulting APs obtained at 1 Hz stimulation frequency. Heterozy- 
gous LQTl mutation prolonged APD by only 8 ms in PFC, in 
contrast to VM, in which a 36 ms APD prolongation is seen. 
Mutation in LQT2 prolongs PFC APD by 28 ms and VM APD by 
46 ms. Clinically, homozygous LQTl mutation leads to the 
Jervell-Lange-Nielsen syndrome[l], which is characterized by 
severe QT prolongation. Homozygous LQTl, simulated as IKs 
deletion, results in further prolongation of VM APD (408 ms 
versus 304 ms WT), which is in stark contrast to the PFC APD 
which exhibits minimal prolongation of 15 ms. 

The differential degree of APD prolongation in PFCs compared 
to VMs is only partially explained by differences in baseline IKr 
and IKs current density. IKr shows a 60% lower whole cell peak 
conductance in PFC, but IKs whole cell peak conductance is 60% 
greater. The lower plateau potential in PFC plays a large role in 
the tissue-specific response via less complete activation of IK 
channels. IKs channel isochronal 1000 ms open probability at the 
PFC plateau potential measures 0.22, compared with 0.82 at the 
VM plateau potential. 

Different potential mechanisms of arrhythmia formation 
in Andersen-Tawil Syndrome 

The Andersen-TawU Syndrome (ATS, or LQTS type 7 [12]) is 
linked to heterozygous loss of function in Kir2.1, which encodes 
the inward rectifier current IKl, through defective channel 
trafficking [13]. Previous work demonstrates that with dominant 
negative suppression of IKl current, delayed afterdepolarizations 
(DADs) are observed in simulated VM[14], and that a PFC focus 
may be responsible for this arrhythmia [15]. 

In Figure 2, dominant negative IKl suppression is simulated in 
both VM and PFC. Arrhythmias occur via different mechanisms 
in the two cell types. In the VM model, progressive downregu- 
lation of IKl leads to DADs of increasing amplitude occurring 
immediately after the AP (panel A), until another AP is triggered 
leading to repetitive membrane depolarization. The mechanism is 
identical to that observed in guinea pig VM[14]. Analogous 
downregulation of IKl in PFC shows no DADs (panel B). The 
slow diastohc depolarization underlying PFC automaticity reaches 
threshold for AP formation earlier than in VM, with a less 
polarized maximum diastolic potential. The result is an approx- 
imate doubling of the intrinsic firing rate stemming from 
unopposed HCN current. Our simulations show that VM carrying 



ATS mutations are capable of triggered activity as an arrhythmic 
mechanism, while corresponding PFC show no DADs but display 
enhanced automaticity, another potential mechanism for arrhyth- 
mia. 

Effect of the SCN5A sodium channel mutations on PFC 
and VM 

SCN5A inactivation defects underlie AP prolongation in LQT3 
via gain of function, producing a persistent component of INa 
termed late current, or INaL. These defects can occur as a result of 
alteration of several distinct biophysical processes. Using several 
different mutations in SCN5A, we explored whether these 
alterations may manifest in a tissue- and/or mutation-specific 
manner. 

The first studied and most thoroughly characterized mutant, 
AKPQ, results in channels that enter a 'bursting' mode of gating: 
channel entry into this mode of gating is rare, but recovery from 
this mode is infrequent, trapping the channel in a state in which it 
cannot inactivate [8] . This results in time-invariant persistent 
channel opening, which generates inward late current that 
prolongs APD. Given the lower plateau potential in PFC, and 
consequent larger driving force for INa entry, our initial 
hypothesis was that a more pronounced effect of AKPQ^ would 
be seen on the PFC APD compared with VM. 

Figure 3 shows the resulting APs using a rate constant for the 
bursting transition that produces 0.5% persistent current, corre- 
sponding to previous studies [8, 16]. At 1 Hz, APD is prolonged by 
69% and 61% in VM and PFC, respectively. After a 2 second 
pause, VM APD carrying the mutation is prolonged to 866 ms 
(increase of 185%), compared to PFC at 892 ms (increase of 
146%). Similarly, more severe APD prolongation is noted in both 
tissue types during bradycardia (pacing at 30 BPM and 40 BPM, 
panels C and D), as summarized in Figure 4E. 

The finding that VM are more affected by AKPQ^ than PFC 
was surprising, since in PFC there is a larger driving force for 
sodium entry during the plateau. We hypothesized that the more 
depolarized diastolic potential in PFC may force a larger fraction 
of SCN5A channels into unavailable states, preventing entry into 
the pathologic bursting states. 

Therefore, we next tested the role of channel availability by 
simulating the LQT3 SCN5A mutation F1473C, a mutation 
characterized clinically by severe QT interval prolongation and 
heavy arrhythmia burden [17], and characterized biophysically by 
a rightward shift in the availability curve of 10 mV and large 
amplitude late current. F1473C channels were simulated with the 




Figure 1. Differential response of VM (panel A) and PF (panel B) to mutation in IKs and IKr channels, leading to LQT1 and LQT2, 
respectively. Compared to wild type (WT, black trace), heterozygous LQTl (blue trace) and LQT2 (red trace), modeled as a 50% reduction in whole 
cell conductance, VIVl APs are prolonged, while PF are mildly affected. Homozygous LQTl (100% reduction in IKs conductance, dashed blue trace) 
leads to substantial APD prolongation in VIVl, and trivial further prolongation in PF. 
doi:1 0.1 371/journal.pone.0097720.g001 



PLOS ONE I www.plosone.org 



2 



June 2014 | Volume 9 | Issue 6 | e97720 



Modeling Tissue Specificity in the Long QT Syndrome 



A Ventricle B Purkinje 




Figure 2. Arrhythmias in Andersen-Tawil syndrome, modeled as dominant negative suppression of IK1 current, have different 
mechanisms in ViVl (panel A) versus PF (panel B) myocytes. For VIVI, compared to WT (blacl< trace) reduction of 1K1 conductance by 90% (blue) 
leads to a large delayed afterdepolarization. Further reduction to 95% (red) increases the size of the DAD, triggering repetitive action potentials. For 
PF, reduction of IK1 conductance maximum diastolic potential is less polarized, leading to increased automaticity. Unstimulated PF APs show basic 
cycle length of automaticity of 2250 msec, which decreases to 1260 with 90% IK1 reduction. 
doi:1 0.1 371 /journal.pone.0097720.g002 



same bursting rate as AKP(]), along with a shift in the availability 
curve of 10 mV in the depolarized direction (see expanded 
methods for values of rate constants) [18]. Figure 4 shows APs at 
1 Hz, demonstrating marked APD prolongation (598 ms in VM 
versus 677 ms in PFC). As with AKPQ, following a 2 sec pause, 
APD is prolonged for both tissue types. Strikingly, for PFC 
following a pause, repetitive early afterdepolarization (EAD) 
formation is observed as a result of severe APD prolongation. 
Similar findings are observed for continuous pacing at slower 
rates. For PFC at 40 BPM, there is a single EAD in every other 
beat (AP alternans, seen in panel D), while at 30 BPM APs with an 
EAD alternate with shorter APs. These findings clearly establish 
an arrhythmogenic substrate within PFC for the F1473C 
mutation. 

S1904L mutation in SCN5A confers a more severe 
phenotype in PFC tlian VM tlirough a different channel 
gating mechanism 

In contrast to typical SCN5A LQT3 mutations, for which the 
phenotype is exacerbated by bradycardia and rest, the SCN5A 
S1904L LCJTS mutation was described in a young patient with 
only a mildly prolonged QT interval who developed palpitations 
with mild exercise [9]. We have recendy characterized these 
mutant channels, which exliibit late current as a result of late 
channel reopenings without a bursting component. 

This channel behavior is recapitulated by reducing the 
inactivation rate constant and reducing the rate constant of entry 
into deeper buried inactivated states by a factor of 5 and 6, 
respectively, generating 0.28% late current (as seen experimental- 
ly[9]). Figure 5A shows the predicted consequence of the S1904L 
in the VM model. As previously shown at 1 Hz [9] , the impact on 
VM is minimal (no prolongation). In contrast, PFC bearing this 
mutation exhibit 120 ms APD prolongation at 1 Hz. More rapid 
pacing (120 BPM) exacerbates this phenotype, prolonging PFC 
APD by 200 ms but stiU with minimal change in VM. The result is 
a markedly flattened restitution curve for the PFC model even with 
a mild increase in rate. A predominant rate-dependent effect on 
PFC explains the emergence of the clinical phenotype with 
exercise as well as the propensity for arrhythmia formation despite 
minimal QT interval prolongation. 



Voltage clamp simulations reveal mechanism of 
mutation-specific differences in LQT3 

The differential response of the PFC AP to AKP(2,and F1473C 
suggests that the shift in channel availability (proarrhythmic) may 
counteract the effect of hyperpolarized diastolic potential in PFC 
(antiarrhythmic). To directly test the contributions of these factors 
individually, voltage clamp simulations were performed at 
diflerent holding potentials (— 100 mV and —80 mV, correspond- 
ing to the diastolic potential of VM and PFC, respectively) with 
steps to different test potentials (+20 and 0 were selected, 
corresponding to plateau potentials). Absolute late sodium current 
at 100 ms was calculated, alongside percent persistent current 
(normahzed to peak) for WT, AKPQ, and F1473C. At a holding 
potential of -100 mV, both AKPQ^ and F1473C show simUar 
amounts of late current, regardless of step potential (and driving 
force), in contrast to WT cells (first two sets in Figure 6A). At a 
holding potential of —80 mV, F1473C exhibits a far greater late 
current than AKPQ. Lesser late current is produced by AKPQ^ at 
depolarized resting potentials due to gating of channels into 
unavailable states (top row of states in Figure IB of Clancy[19]) 
from which they cannot burst nor conduct current during the 
plateau of the AP. Thus, the higher driving force for current entry 
is counterbalanced by lower fractional occupancy of bursting states 
that carry INaL. When availability is shifted in F1473C channels, 
direct comparison of the -100 to +20 mV simulation (mimicking a 
VM AP profile) and —80 to 0 mV simulation (mimicking a PFC 
AP profile) shows more absolute late current in the PFC-like 
profile. This substantial increase in late current underlies the more 
severe prolongation of APD in PFCs for the F1473C mutation. 

We also used the square voltage clamp protocol for the S1904L 
mutation to mimic VM and PFC AP morphologies. In the "VM- 
like" voltage clamp, the magnitude of the late current (and late 
current fraction) is small compared to the PFC-like voltage clamp 
profile (see Figure 6B). The prolonged time course of inactivation 
at 0 mV compared to +20 mV (time constant of current decay 
67 ms at 0 mV versus 30 ms at +20 mV) confers a larger 
functional impact in the context of the S 1 904L mutation, allowing 
persistent channel opening during the step potential. This 
corresponds to 85% smaller net summed current entry (net 
charge) during the test pulse for a VM-like profile (bar graph. 
Figure 6B). Repetition of these voltage clamp waveforms at 2 Hz 
does not appreciably alter net current entry, suggesting that the 
larger effect at higher frequencies in PFC is due to smaller 
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Figure 3. VM and PFC carrying tKie SCN5A AKPQ mutation exhibit similarly severe APD prolongation. Panel A and B) Effect of a 2 
second pause for ViVI and PFC, respectively, compared with 1000 msec (last beat of drive train shown for WT, black trace, and mutant cells, blue 
trace). Panels C and D) Steady state rate adaptation for VIVI and PFC, respectively (30 BPM: blue trace, 40 BPM: red trace, 60 BPM: black trace). E) 
Summary results for APD prolongation for each cell type for different stimulus protocols; VM: black bars, PF: gray bars. 
doi:1 0.1 371 /journal.pone.0097720.g003 



counterbalancing net repolarizing current such as the transient 
outward current. 



Discussion 

The major findings of this study are that a) LQTl and LQT2 
exert a relatively mild phenotype in PFC; b) different mechanisms 
of arrhythmia formation are seen in Andersen-Tawil syndrome 
based on cardiac tissue type; c) the canonical LQT3 mutant 
AKPQhas severe phenotype in both VM and PFC tissue types; d) 



in contrast to AKPQ, the LQT3 mutant F1473C shows a 
particularly severe phenotype in PFC with repetitive membrane 
depolarizations due to shifted availability of sodium channels; and 
e) the LQT3 mutant S1904L exerts a mild phenotype, but is more 
severe in PFC with strong rate-dependent exacerbation of 
phenotype. Numerous lines of evidence implicate the PF system 
as important in the genesis of LQTS-related arrhythmias, 
including findings from electrophysiological studies of affected 
patients. Our findings, which stem from inherent differences in the 



PLOS ONE I www.plosone.org 



4 



June 2014 | Volume 9 | Issue 6 | e97720 



Modeling Tissue Specificity in the Long QT Syndrome 



Ventricle 




B 



Purkinje 



WT 

F1473C 



L 



500 ms 




F1473C rate adaptation 

D 




60BPM 
40BPIV1 
30BPIVI 



250 ms 




Figure 4. Differential effect on APD for VWl versus PF of tlie SCN5A F1473C mutation. Panels A and B) Effect of a 2 second pause for VIVI 
and PFC, respectively, compared with 1000 msec (last beat of drive train shown for WT, black trace, and mutant cells, blue trace). APD in the F1473C 
shows pause-dependent increase, with spontaneous oscillations in membrane potential from sequential early afterdepolarizations seen post-pause 
for PF but not VIVI. Panels C and D) Rate adaptation for VM and PFC, respectively (30 BPM: blue trace, 40 BPM: red trace, 60 BPIVI: black trace. In VIVI 
cells, substantial APD prolongation is seen (on the order of that seen with AKPQ mutation) with pronounced bradycardia dependence. In PF APs, 
severe APD prolongation is seen at 60 BPM, for 40 BPM stable alternans are seen (with a small early afterdepolarization on every other beat), and at 30 
BPM the membrane fails to repolarize. 
doi:1 0.1 371/journal.pone.0097720.g004 



electrophysiology of the PFC compared to the VM, provide a 
theoretical basis with which to understand these chnical observa- 
tions. 

Arrhythmias in the ATS may result from enhanced PFC 
automaticity 

The mechanism of arrhythmia formation in ATS has not been 
estabhshed. Previous simulation studies showed that in VM with 
dominant negative suppression of IKl, DADs were ob- 
served[14,20]. Our findings in VM reproduce this behavior, with 
repetitive membrane depolarization upon critical IKl reduction 
(Figure 3). PFC with dominant negative IKl suppression, 
however, did not show evidence of DADs. Instead, we saw a 
graded increase in intrinsic automaticity. This suggests that 
abnormal automaticity (stemming from the cardiac conduction 
system) and/or triggered activity (from VM) may underlie the 
ventricular arrhythmias seen in ATS. The Sanders lab has 
reported sharp "Purkinje potentials" preceding ventricular 
premature beats during electrophysiology studies in ATS pa- 
tients [15]. The bidirectional ventricular tachycardia often seen in 
these patients may reflect alternation between two different PFC 
foci with enhanced automaticity, via a mechanism similar to that 
described by Weiss [21]. Our results thus provide a plausible 
mechanism that may underlie the chnical demonstration of His- 
Purkinje system involvement. 



The importance of the plateau potential in determining 
APD for LQT1 and LQT2 

In all mutations studied, one common theme is the importance 
of the plateau potential (and plateau current magnitudes) in 
determining the effect of a mutation. Our results simulating LQT 1 
and LQT2 (Figure 1) show a more profound APD-prolonging 
effect of IK deletion in VM. It has been widely reported that 
baseline IK density is different in VM and PFC [22], and as such, 
the dependence on IK for membrane repolarization (and 
consequence of deletion of these currents) would be expected to 
be different in each cell type. However, we found that differences 
in channel activation may underlie the greater APD prolongation 
in VM. Indeed, at the VM plateau potential, our results predict a 
nearly four-fold greater IKs channel open probability. Previous 
experiments using APs as command waveforms predicted a similar 
result, wherein larger IK current was elicited using an AP with 
ventricular morphology [2 3]. Taken together, these results suggest 
that not only channel density, but AP shape regulates the effect of 
LQTl and LQT2 mutations in different tissue types. 

Our model for LQTl was relatively simple, in that a simple 
change in conductance was simulated rather than a more subtle 
loss-of-function resulting from altered kinetics. Pourrier and Nattel 
reported that accessory subunits of the IKs channel may express 
differentially in different cardiac tissue types, with a higher density 
in Purkinje fibers [24]. A disproportionate reduction in PFC IKs 
via differential expression of subunits may modulate the predic- 
tions of the simpler model described herein. 
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Figure 5. The S1904L mutation, characterized clinically by 
arrhythmia at low levels of exertion, simulated at clinically 
relevant pacing frequencies. The mutation exhibits minimal 
response in VM (panel A, traces overlap) and greater APD prolongation 
in PF (panel B), with APD alternans seen at 120 BPM (lower row). WT 
cells: black trace, S1904L: blue trace. 
doi:1 0.1 371 /journal.pone.0097720.g005 

Altered Na+ channel availability and defects in 
inactivation modulate tissue-specific effects of LQT3 
mutations 

Clancy [8] showed that AP prolongation with the AKPQ^ 
deletion in SCN5A can arise from sodium channels that enter 
an abnormal bursting mode, in which the channel is trapped in a 
configuration fi-om which it cannot transition to inactivate. 

We initially hypothesized that for PFC carrying AKPQ^ a larger 
driving force for sodium current entry during the relatively 
hyperpolarized plateau provides the basis for more severe APD 
prolongation. We did not find this to be the case, instead observing 
a slightly more pronounced APD prolongation in VM compared 
to PFC. The steady-state availability curve for SCN5A predicts 
0.68 fractional availability at the PFC diastolic potential of — 
83 mV and 0.85 fractional availability at the corresponding VM 
potential of —91 mV. The additional fraction of inactivated 
channels in PFC cannot transition into the bursting states to 
produce persistent current. 

Our previous work described the LQT3 mutation F1473C, 
characterized by severe QT prolongation that is further exacer- 
bated by bradycardia and substantial late sodium current [18]. 
When implemented into the PFC model, the results indicate that 
the shift in availability seen with F1473C compensates for the 
relatively depolarized resting potential seen in this tissue type. 
With the consequent greater availability of mutant channels, the 
enhanced driving force for sodium current entry in PFC causes 
disproportionate prolongation of APD. Moreover, we demon- 
strated in this study that PFC are prone to EADs, the putative 



cellular triggers of ventricular arrhythmia, and stable AP alternans 
promoted by slow pacing rates. These findings suggest that 
disruption of PFC electrophysiology was responsible for the 
affected patient's unusually heavy burden of ventricular arrhyth- 
mia clinically, which required defibrillator implantation as a 
neonate [17] and which improved markedly upon faster pac- 
ing[25]. 

Our finding that channel steady-state availability impacts 
channel function in a tissue-specific manner led us to explore 
whether our results have important implications for drug therapy. 
We tested two different local anesthetic drugs, flecainide (modeled 
as a charged local anesthetic SCN5A blocker) and lidocaine 
(modeled as an uncharged local anesthetic SCN5A blocker which 
affects channel availability) on our PFC model of F1473C. In the 
drug free state, severe APD prolongation is seen along with 
multiple EADs, particularly in response to a 2 second pause. 
Flecainide was modeled as a 20% block in sodium current and 
lidocaine as a 20% block in sodium current plus a leftward shift of 
the availability curve of SCN5A. We identified that with lidocaine, 
the availability shift of lidocaine conferred greater efficacy in terms 
of APD shortening and abolishing EADs (unlike flecainide). These 
simulations raise the intriguing possibility that drug efficacy may 
vary in a tissue-specific manner in LQTS patients. 

The LQT3 mutation S1904L is unique in that APD prolon- 
gation and ventricular arrhythmias occur with faster heart rates 
(unlike in AKPQ,or F1473C). Furthermore, arrhythmias occurred 
with mild exercise in the proband carrying this mutation [9]. Our 
previous characterization of this mutation showed minimal effect 
on ventricular APD at physiological pacing frequencies (although 
an effect was measurable at >200BPM[9]). In our current study, 
we show that VM exhibit minimal change in AP duration at both 
60BPM and 1 20BPM, whereas for PFC a more pronounced effect 
is seen. Via this mechanism, a heterogeneous repolarization 
substrate emerges preferentially at higher heart rates, explaining 
the propensity toward arrhythmia formation during exercise. 

Impact of concomitant repolarizing current blockade on 
LQTS phenotype 

Given its relatively slow activation, IK in cardiac cells is thought 
to play a larger role under conditions that cause AP prolongation. 
Acquired LQTS, as occurs via partial blockade of IKr associated 
with QT-prolonging drug therapy, may reflect unmasking of a 
subclinical channelopathy, through deletion of an important 
repolarizing influence that in WT cells may be redundant. As 
studied by Gintant in canine Purkinje fibers [26], AP prolongation 
from repolarizing current blockade may manifest prominently in 
PF. Since there are known difficulties in estimating the values of 
IK density in VM, and this density may vary from ceU to ceU, we 
explored whether variation in IK in cells carrying LQT3 
mutations may affect our findings. As shown in Figure S3 in File 
SI, we found that concomitant IK block leads to substantial AP 
prolongation in VM carrying LQT3 mutations (both AKPQ_ and 
F1473C), but no EADs are seen even with IK block up to 50%. 
This suggests that while repolarizing currents play a larger role in 
the presence of coincident gain-of-fuiictioii LQTS mutations, it is 
the degree of late current augmentation that determines whether 
triggered activity occurs. 

We furthermore performed simulations in another model of the 
human ventricular action potential, and with both epicardial and 
endocardial profiles, to establish that our findings in LQT3 
generalize across tissue layer and do not reflect artifactual behavior 
unique to one ionic model of the ventricular myocyte. For these 
simulations (described in further detail in File SI), the model of 
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Figure 6. Voltage clamp protocols demonstrate the mechanism of late current generation in LQT3 mutations. Panel A illustrates the 
voltage clamp waveform used. Panel B demonstrates the important role of resting membrane potential and plateau potential on degree of late 
sodium current in AKPQ (red bar) and F1473C mutations (blue bar) compared to WT (black bar). For a "VM-like" voltage clamp protocol (holding - 
100 mV, step potential 20 mV, top row), similar persistent current is seen for KPQ and F1473C mutations. When holding potential is changed to — 
80 mV mimicking PF AP resting potential, second row, persistent current is less for the KPQ mutation than for the F1473C mutation (which has 
increased availability at this potential). Third and fourth rows show the individual effect of altering plateau and resting potentials. Panel C) Response 
of step potentials simulating VIVl (black) and PF (blue) in cells carrying the S1904L mutation. The "PF-like" voltage clamp waveform produces more 
inward current at all time points, leading to more accumulated inward charge as demonstrated below in the bar graph. 
doi:1 0.1 371 /journal.pone.0097720.g006 



Ten Tusscher[27] was utilized. Tlie fmding.s confirm those 
obtained using the model described here (Figure S4 in FUe SI). 

Nature of the interaction between PFC and VM impacts 
arrhythmogenesis at the organ level and potential 
therapeutic strategies 

The cell type-specific responses to different LQTS mutations 
theoretically provide a substrate for arrhythmogenesis, by increas- 
ing the heterogeneity of repolarization (contributing to dispersion 
of APD and dispersion of the QT interval). Evidence is 
accumulating that the PF network is organized in complex sheets 
adjacent to the endocardial surface of the heart that only interact 
with VM at sparse PFC-VM junctions. This geometry may 
facilitate propagation of EADs that occur preferentially and focally 
in PFC to neighboring cell layers, suggesting the potential for 
reentry and proarrhythmia. 

Antiarrhythmic therapy is the standard of care in the 
management of many patients with LQTS, with an adjunctive 
role of implantable cardioverter-defibrillators to reduce mortality 



associated with sudden cardiac death. A central role of the PFC 
and specialized conduction system may suggest specific pharma- 
cologic targets that may be tailored to a specific mutation. 
Furthermore, since PFC can be targeted anatomically, alternate 
modes of therapy, for example, via catheter ablation may be 
possible. We have recendy characterized a PFC-specific mutation 
in the transient outward current, identified in patients with 
idiopathic ventricular fibrillation who were successfully treated 
with catheter ablation [28]. Since the His-Purkinje system inserts 
endocardially, PF triggers in other heritable channelopathies may 
evolve to represent attractive targets of ablation therapy, as has 
recently been reported in a canine model of LQT3[29] and during 
electrophysiology studies of patients with LQTS [2]. 

Limitations 

A first limitation of the simulations shown is that conclusions 
apply to the single cell, but without multicellular analysis, its 
impact on the whole organ level remains untested. Future work by 
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our group will explore the potential for triggering of arrhythmias 
from these cellular proarrhythmic states. 

A second limitation lies in the formulation of calcium handling 
utilized. This Ls a weakness of ever)' current model of the PFG. We 
plan to update our model to incorporate what is known about PFC 
ultrastructurc and more accurately reflect intracellular calcium 
cycling. Such an addition may further elucidate the link between 
altered membrane repolarization and triggers of arrhythmia. 

Conclusions 

The Purkinje fiber network has long been implicated as a source 
of ventricular arrhythmia. We show that for several mutations 
linked to LQTS, a different phenotype emerges in difierent 
cardiac tissue types, often predicting a more severe phenotype in 
PFC. Our findings suggest the potential for a new type of 
regulation of heritable channelopathies in which not only channel 
dysfunction but the variable impact on different cardiac tissues can 
modulate the severity of phenotype, with important implications 
for therapy of these conditions. 

Methods 

The whole-cell mathematical models utilized in this study have 
been described previously [1 0, 1 1]. Briefly, these models are 
composed of systems of coupled nonlinear diff("r("ntial equations 
representing channel occupancies, which are integrated numeri- 
cally using a fourth order Runge-Kutta ordinary differential 
equations solver with an adaptive time-step. Computations were 
performed in C-H- language on a 4-processor workstation, 
allowing single cell computations in real time. Further details 
regarding the computational models can be found in previous 
work[10,ll]. 

The previously developed PFC model[l 1] is implemented 
without change. In order to model the mutations identically in 
each tissue type, the Markov representations of membrane 
currents from the VM model were replaced with the correspond- 
ing models from the PFC, specifically the slow component of the 
delayed rectifier potassium current (IKs) and the fast sodium 
current INa[l 1]. To ensure appropriate balance of densities of the 
fast component of the delayed rectifier IK (IKr) and IKs after this 
change, conductance of these currents was adjusted to produce 
appropriate response of APD to specific blockers of each 
component (E-4031 for IKr[30] and chromanol 293B for 
IKs [31]). Appropriate APD for human epicardial VM was 
reproduced, chiefly through a 20% increase in the conductance 
of the L-type calcium current (matching recent data recording L- 
type current density in human right[32] and left[33-35] VM). A 
complete description of model parameters is provided in the 
expanded methods section in the File SI. The model reproduces 
key features of human epicardial VM, including appropriate rate 
dependent shortening of APD, as seen in Figure SI in File SI. In 
Figure S2 in File SI the updated model output is compared to a 
recorded action potential and the original model of Iyer et al [10], 
showing minimal difference in action potential shape and 
duration. 
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